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Abstract 

The vibration ring was conceived as a driveline damping device to prevent structure-borne noise in 
machines. It has the appearance of a metal ring, and can be installed between any two driveline 
components like an ordinary mechanical spacer. Damping is achieved using a ring-shaped piezoelectric 
stack that is poled in the axial direction and connected to an electrical shunt circuit. Surrounding the stack 
is a metal structure, called the compression cage, which squeezes the stack along its poled axis when 
excited by radial driveline forces. The stack in turn generates electrical energy, which is either dissipated 
or harvested using the shunt circuit. Removing energy from the system creates a net damping effect. The 
vibration ring is much stiffer than traditional damping devices, which allows it to be used in a driveline 
without disrupting normal operation. In phase 1 of this NASA Seedling Fund project, a combination of 
design and analysis was used to examine the feasibility of this concept. Several designs were evaluated 
using solid modeling, finite element analysis, and by creating prototype hardware. Then an analytical 
model representing the coupled electromechanical response was formulated in closed form. The model 
was exercised parametrically to examine the stiffness and loss factor spectra of the vibration ring, as well 
as simulate its damping effect in the context of a simplified driveline model. The results of this work 
showed that this is a viable mechanism for driveline damping, and provided several lessons for continued 
development. 


Nomenclature 


A Piezo element face area 

c Dynamic compliance of system 

C Piezo stack capacitance with free mechanical boundaries 

d Piezo constant of element 

^ e ff Effective piezo constant of stack 

D Electric displacement of piezo element 

E Electric field of piezo element 

/ Frequency 

f n Natural frequency of system with vibration ring 

fn,b Natural frequency of baseline system 

F Piezo stack force 

F Piezo stack force with reversed sign convention (compressive force > 0) 
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Compression cage force 

Piezo element force 

Vibration ring radial force 

System input force, applied to shaft 

Vibration ring link force 

System output force, transmitted to housing 

Force transmissibility of system 

Force transfer ratio, efficiency of force transfer from vibration ring to stack 

Function relating complex stack stiffness to complex vibration ring stiffness 

Function relating stack short circuit stiffness to complex stack stiffness 

The value of H at the system natural frequency 

Piezo stack current 

Piezo element current 

Load current 

Imaginary number 

Piezo stack stiffness 

Piezo stack complex stiffness 

Baseline system stiffness 

Baseline system complex stiffness 

Compression cage stiffness 

Effective complex stiffness of system (with vibration ring) 

Effective stiffness of system at the natural frequency 

Piezo stack stiffness at the system natural frequency 

Open circuit stiffness of piezo stack 

Short circuit stiffness of piezo stack 

Vibration ring stiffness 

Complex valued vibration ring stiffness 

Vibration ring stiffness at the system natural frequency 

Piezo element thickness 

Length of mechanical amplifier links in vibration ring 
Axial length of vibration ring hoop walls 
Mass of shaft 

Number of piezo elements in a stack 
Power delivered to load 
Load resistance 

Load resistance associated with maximum piezo stack loss factor 
Elastic compliance constant of piezo element under constant electric field 
Piezo element strain 
Time 
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Piezo element stress 
Piezo stack voltage 
Piezo element voltage 

Distance between inner and outer vibration ring links 

Electric energy delivered to load during electromechanical coupling virtual 

experiment 

Mechanical energy recovered during electromechanical coupling virtual 
experiment 

Energy dissipated per steady state oscillation cycle of a complex spring 
Total steady state oscillation energy of a complex spring 
Radial compression of vibration ring wall 

Unloaded radial distance between inner and outer walls of vibration ring 
Loaded radial distance between inner and outer walls of vibration ring 
Shaft displacement 
Piezo stack displacement 

Piezo stack displacement with reversed sign convention (compression > 0) 
Piezo element displacement 
Unloaded height of piezo stack 
Loaded height of piezo stack 

Maximum stack displacement in the electromechanical coupling factor virtual 

experiment 

Load impedance 


a Ratio of cage stiffness to vibration ring stiffness 

P Ratio of system natural frequencies with and without vibration ring 

e t Permittivity constant of piezo element under constant stress 

h Piezo stack loss factor 

r| max Maximum loss factor of piezo stack 

r\b Loss factor of baseline system 

^Ivr Vibration ring loss factor 

r|^. ax Maximum loss factor of vibration ring 

k Piezo stack electromechanical coupling factor 

k c i Piezo element electromechanical coupling factor 

^unload Unloaded vibration ring link angle 

<|>load Loaded vibration ring link angle 

co Angular frequency 
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1.0 Introduction 


Rotating machines create noise with levels ranging from unpleasant to damaging. The main noise 
sources are high speed driveline components, such as gears, bearings, shafts, fans, etc. These parts create 
periodic forces, which shake the machine body and generate structure -borne noise. The example that 
motivated this work was the environment within the cabin of a rotorcraft. Here the driveline-induced 
noise exceeds 100 dB, making it impossible to communicate without headsets. This environment creates 
pilot fatigue and limits the commercial use of rotorcraft for civilian transportation. As a part of the 
National Aeronautics Research and Development Plan, the United States set the goal to reduce main rotor 
gearbox noise by 20 dB (National Science and Technology Council, 2010). 

Structural damping is the most common method to reduce machine noise, accomplished by applying 
hysteretic material (like rubber) to the radiating surfaces. As the surfaces move, work is done on the 
material, which removes part of the energy in the form of heat. For rotorcraft this approach is limited, 
because it would require too much added mass to effectively dampen the cabin walls. This shortcoming 
motivates the present investigation to develop a damper that would reduce vibration energy within the 
driveline. The idea is conceptually illustrated in Figure 1 using the case of a simple gearbox. Here 
dampers have been placed around the shaft bearings to reduce the amount of vibration energy transferred 
to the gearbox housing and to the rest of the machine. 



Figure 1 . — Concept drawing of driveline dampers mounted around the shaft bearings of a 
gearbox. The dampers reduce the driveline vibration and in-turn lower the amount of 
energy transferred through the gearbox housing to the rest of the machine. 
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Figure 2. — Simple driveline with and without in-line damper, (a) Schematics, (b) Analytical models. 


The damper was envisioned to have a simple ring-shaped geometry, allowing it to be placed in- 
between any two driveline components (e.g., between a gear/shaft or bearing/housing). In-line damping 
devices are not typical, because they alter the frequency response of the original system. To gain some 
intuition for how to best configure such a device, consider the simplified driveline schematic in 
Figure 2(a). On the left is a shaft supported by a bearing within a mechanical housing. On the right, is the 
same system including the subject in-line damper. Also depicted are analytical models, in Figure 2(b). 
Here the shaft is represented by a mass and the bearing and damper are represented by springs with 
stiffness (K) and loss factor (r|) parameters, define as follows. 

k b = k b (l + jr\ b ) (p 

Kx =Kx (l + /nvr) (2 

The series combination of springs is expressed as, 


^eff - “ 




k h + K 


(3) 


To compare these two driveline cases, the metric of force transmissibility (FT), 

pT _ ^out ^eff 

F m ~(2nff m s +k eS 


(4) 


is used to express the transfer of force from the shaft to the housing in the frequency domain. The results 
are plotted in Figure 3. With no damper, there is a sharp peak in the FT at the natural frequency. As 
shown in Figure 3(a), the peak is reduced by lowering the damper stiffness; however, this simultaneously 
reduces the natural frequency. To avoid drastic frequency shifts, the damper must be significantly stiffer 
than the bearing. Figure 3(b) shows the case where the damper stiffness is fixed at two times the bearing 
stiffness and the loss factor is varied. Here, it’s evident that increasing the loss factor provides the desired 
damping effect without substantial frequency-shift. 
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Figure 3. — Force transmissibility plots showing effects of an in-line damper, (a) Effect 
of varying damper stiffness with a constant damper loss factor, (b) Effect of varying 
damper loss factor with a constant damper stiffness. 


Traditional damping materials, such as rubber, exhibit high loss factor but are far too soft to be placed 
in-series with steel driveline components. Shunted piezoelectric (piezo) stacks (Hagood & von Flotow, 
1991) were found to be a suitable alternative. Piezo stacks are essentially multilayer piezo devices, which 
generate significant electric charge in response to applied vibratory force. By shunting their electric 
terminals the generated electrical energy is removed from the system to create a net damping effect. 

When used in compression they have an elastic modulus comparable to steel with loss factor comparable 
to rubber. There are also a few practical limitations associated with the technology. The damping effect is 
primarily limited to the axis in which the stacks were electrically poled. Piezo stacks have high 
compressive strength along their poling axis, but are relatively fragile when stressed in tension or shear. 

A survey of the literature revealed that shunted piezo stacks have been experimentally evaluated for 
in-line driveline damping, as reported in (Atzrodt, Mayer, & Melz, 2009). The test setup consisted of a 
single shaft supported by a bearing, as illustrated in Figure 4. Four piezo stacks were placed in-between 
the bearing and mechanical housing, and each was electrically shunted by a resonant circuit. Radial 
vibratory force was applied to the shaft using a shaker, and the resulting shaft velocity response spectrum 
was examined. Engaging the resonant circuits reduced the response at the natural frequency by 17.5 dB. 
This result provides strong support for the use of shunted piezo stacks. However, this particular stack 
arrangement appears to be limited to use between a bearing and housing. On a rotating part the loads 
would be distributed, requiring the use of numerous stacks to create an even reaction. In addition, torque 
would be transferred through the stacks creating the potential for shear failure. 
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Piezo stacks 


Figure 4. — Depiction of piezo shunt damping 
experiment reported in (Atzrodt, Mayer, & 
Melz, 2009). 



Figure 5. — Original concept drawing of the vibration ring. 


The vibration ring was conceived as a mechanism to implement piezo shunt damping anywhere 
within the driveline. As depicted in Figure 5, it’s shaped like a ring and could be installed between any 
two driveline components like an ordinary mechanical spacer. Instead of using radially oriented piezo 
stacks, a single ring-shaped stack is used. Radial force applied to the mechanism at any position is 
transferred to the stack through the compression cage. This generates charge in the stack, which is 
dissipated by the electric circuit to create the damping effect. The compression cage is the key to this 
innovation. It has a unique concave structure that squeezes the piezo stack in the axial direction when the 
inner and outer walls of the vibration ring are moved closer together. The compression cage also holds the 
stack under a compressive preload, and effectively shunts torque, to protect the stack from potentially 
damaging tensile and shear forces. 

This new vibration control mechanism is being evaluated in a two-phase project supported by 
NASA’s Aeronautics Research Mission Directorate Seedling Fund. This report provides a comprehensive 
record of phase 1 , where the feasibility of vibration ring was examined through a combination of design 
and analysis. 
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2.0 Design 

Preliminary design work focused on finding a suitable compression cage structure. Several iterations 
of solid modeling and rapid prototyping were done to find a concept that effectively converted radial 
motion into axial motion, and could also be manufactured and assembled. Once a basic understanding of 
the design challenges was reached, rough guidelines were formulated to provide more structure to the 
design process. Subsequent design concepts were developed with more detail, including finite element 
analysis (FEA) and complete prototype development. Test hardware was also created to evaluate the 
damping properties and power generation of the prototypes. 

2.1 Preliminary Design Concepts 

Though there were several preliminary designs created, only select concepts are discussed here to 
emphasize lessons learned or to provoke further consideration. As a first step in the design process, 
several piezo device manufacturers were contacted to find a ring-shaped piezo stack. Most offered to 
manufacturer custom designs, but none had this stack geometry available as an off-the-shelf product. One 
company provided a set of piezoceramic rings, shown in Figure 6(a), as building blocks to create a stack 
in-house. A second company suggested combining several commercially available linear stacks using 
metal rings or “stack rings”, such as in Figure 6(b). The preliminary vibration ring design concepts were 
based on the geometry of the piezoceramic rings. The rings were found to be too small to be practical, so 
the detailed designs discussed in the subsequent section made use of the stack-ring concept. This is 
considered a temporary solution until a suitable ring-shaped stack can be procured. 

Concept 1 (Coni) was created to mimic the original depiction of the vibration ring. As shown in 
Figure 7, a solid model of the compression cage and corresponding plastic prototype were made. The 
cross section of the compression cage looked like a set of wings, and so this profile was given the name 
‘butterfly amplifier’. By squeezing the walls of the plastic prototype at the open end the axial 
compression action of the butterfly amplifier could be observed. One of the major issues identified was 
that the design required significant axial width to accommodate the stack and the butterfly wings. To 
illustrate this, the model and prototype are shown in the figure with two bearings stacked axially. 



Figure 6. — Stack options, (a) Piezoceramic rings, (b) Stack ring concept. 



Figure 7. — Coni, (a) Solid model, (b) Plastic prototype. 
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Subsequent attempts were made to redesign the compression cage, so that it could be assembled 
around the stack. However, none of the ideas evaluated were practical. In Con6 a temporary solution was 
implemented so that other aspects of the design could start to be explored. The solid model is shown in 
Figure 8(a) and corresponding plastic prototype in Figure 8(b). Here, the compression cage was made in 
two axial pieces that were assembled around the stack using small bolts. Unlike Coni, the plastic cage 
parts felt rigid and the compression action could not be observed. The difference was that Coni was cut 
along its cross-section allowing for more freedom of movement at the boundary. This led to the 
realization that the compression cage needed to be segmented circumferentially to function correctly. 

It was felt that additive manufacturing would have some advantages for developing metal prototypes 
of the compression cage. The prototypes could be made directly from solid models (without drawings). 
Also, the compression cage could be built with complex 3D geometry that would otherwise be impractical 
to machine. To begin exploring this approach, the Con6 compression cage was built out of stainless steel 
using the direct metal laser sintering (DMLS) process, as well as by investment casting. One piece of each 
cage structure is shown in Figure 9. The DMLS part was found to be more precise than the casting. 
However, the casting had smoother surfaces and less porosity. The relative porosity of the parts was 
inferred, as the DMLS part was lighter and had a lower frequency acoustic response when struck. 



Figure 8. — Con6. (a) Solid model, (b) Plastic prototype. 


(a) Sintered (b) Casted 



Figure 9. — Con6 stainless steel parts, (a) Sintered, (b) Casted. 
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Con8 and 9 did not contribute directly to the development of the more advanced designs, but they are 
discussed here to highlight specific features. Con8, illustrated in Figure 10, was a first attempt at creating 
a segmented cage. For each segment, the inner and outer wings were separate parts. Theoretically the 
compression cage could be assembled around the ring-shaped stack using a simple three-step process. 
First, two metal rings would be placed above and below the stack. Second, each wing segment would be 
inserted radially into slots in the metal rings. Third, the inner and outer hoops would be slid onto the 
design axially. 

Con9, as shown in Figure 11, used an expanding form of the compression cage structure that was 
totally different than the other concepts. The first version, Con9a, is shown in Figure 1 1(a). The green 
section provided the amplifier action. Both the external and internal vertical members have a c-shape to 
them, but with different degrees of curvature. Radial deformation of the external member is transferred to 
the internal member through a rigid link. The internal member in turn straightens and compresses the 
stacks that are above and below it. Con9b, shown in Figure 1 1(b), was a symmetric version of this design 
that also had flat external surfaces for mating with other parts. The expanding structure of Con9 was 
interesting, but it required more axial width than the butterfly amplifier, it did not have an assembly 
solution, and it did not appear to offer any performance benefits to justify further investigation. 



Figure 11. — Con9 solid models, (a) Con9a, (b) Con9b. 
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TABLE L— DESIGN SPECIFICATIONS 


External dimensions 

ID = 52 mm 

OD < 86 mm 

Axial width 15 to 30 mm 

Load 

Radial compression 
Nominal 400 N 
Maximum 845 N 

Stack type 

Noliac piezoceramic stack NAC2013-H16-A01 
Height 15.8 mm 
Face 5.0 by 5.0 mm 

Static FEA performance criteria 

Maximize the force transfer ratio at nominal load 



Photograph of stack 


2.2 Detailed Designs 

After evaluating the preliminary concepts, the lessons-leamed were applied to create more detailed 
and functional designs. The new designs used the butterfly amplifier with a segmented configuration. 

Also, the stack-ring concept from Figure 6(b) was implemented as a temporary substitution for a ring- 
shaped stack. Specifications for the external dimensions, load, stack type, and FEA performance criteria 
were established, as listed in Table 1. In anticipation of future system level testing, the dimensions and 
load specifications were made to correspond to specific driveline test rigs at NASA. For this preliminary 
work, only compressive radial loads were considered. The stack type and FEA criteria were chosen using 
guidance from the analytical framework. For damping at the stack level, a high electromechanical 
coupling factor was shown to be essential (Eq. (31)). The specified stack was selected because it had a 
relatively high value (0.72) amongst commercially available products. To transfer the stack damping 
properties to the vibration ring, the framework showed that the ratio of cage stiffness to vibration ring 
stiffness needed to be minimized (Eq. (50)). As a means to achieve this goal indirectly, the FEA 
performance specification called for maximizing the force transfer ratio (Eq. (53)) between the vibration 
ring and stack. It should be noted that only static FEA was implemented for this tuning process. 

ConlO is illustrated in Figures 12 to 14. This design had few parts, but required additive 
manufacturing to be created. Static FEA was used to evaluate the loaded response, in terms of stiffness, 
force transfer ratio, and stress. This is described in detail in the analysis section of this report. The first 
version, ConlOa, is shown in Figure 12. It used nine stacks and had four sections: Top and bottom pieces 
to the cage structure, as well as inner and outer hoops. To assemble this design, the bottom cage section 
would be bolted to the inner and outer hoops, the stacks would be inserted and then the top section would 
be bolted into position. FEA results showed that the maximum stress was highly dependent on the link 
radii and thickness, and that the design would be overly sensitive to fabrication error. This was attributed 
to the tight radius in the butterfly wing, as annotated in Figure 12(b). In addition, the bolt clearances were 
considered to reduce assembly precision. 

Con 10b was a modified version, intended to overcome the limitations. It was a two piece structure 
consisting of the main compression cage section, shown in Figure 13(a), and an outer hoop (not shown). 
The piezo stacks were intended to be inserted radially, the hoop would be slid around the structure, and 
then countersunk screws would be used to precisely fasten the two members. Relative to ConlOa, the 
fasteners were located further outboard to eliminate the tight radius in the butterfly wing. This design was 
still somewhat sensitive to link geometry, and so the DMLS process was selected over investment casting 
to create the prototype. 

All of the thin stainless steel sections for Coni Ob were severely warped by heat in the DMLS process. 
To allow the parts to mate together, quite a bit of final machining and bending of the parts needed to be 
done. Shims were then inserted to apply preload to the stacks. The final prototype is illustrated in 
Figure 14. 
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Figure 12. — ConlOa solid model, (a) Overview, (b) Cross section. 



Figure 13. — ConlOb solid model, (a) Main compression cage section, (b) Cross section. 



Figure 14. — ConlOb prototype created by DMLS. Shown with an 
aluminum mounting fixture inside. 
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Coni 1, shown in Figures 15 to 18, was developed to improve the prototype precision by using more 
conventional manufacturing processes. The first iteration, Coni la, is illustrated in Figure 15. Instead of 
having a unified cage structure, there were several separate parts that could be made using a combination 
of wire electrical discharge machining, milling, and turning. The parts installed axially, and were loaded 
against metal tabs on each side. A couple of additional special features were added to make assembly 
easier. First, the links were notched in the center, so they could be folded slightly and inserted between 
the inner and outer hoops. Once in place, the links would be released and would expand into small 
pockets in the hoop walls. The pockets were essentially a locating feature. The stack cap part would fill 
the notch in the link, like a tongue-in-groove arrangement. Second, set screws were included in the tabs to 
hold the assembly in place, and apply preload to the stacks. 

The pocket and tongue and groove features of Coni la were thought to complicate fabrication, and so 
Coni lb was created without these features. As illustrated in Figure 16, the design was otherwise very 
similar to Coni la. One notable change was that the stack rings were placed outside of the links to allow 
them to be removed. This provided the option to measure their influence on performance. 




* 



t c i M 




Stack 

ring 

Pocket 


Figure 15. — Coni la solid model, (a) Exploded view, (b) Cross section. 



Figure 16. — Coni 1b solid model, (a) Overview, (b) Cross section. 
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Manufacturing of the Coni lb prototype was fairly straightforward. There were no issues with 
tolerances or warped parts. However, assembly was more difficult than anticipated. Two pictures from the 
assembly process are included in Figure 17, for reference. The main issue was that the interference fit 
between the links and hoops made it so the parts needed to be forced into position. To install the final 
links, the hoop walls needed to be pried apart, as illustrated in Figure 17(b). After all of the parts were 
installed, the inner and outer hoops were slightly offset in the axial direction. This occurred, because the 
links had rotated to relieve their interference fit. To correct this, a mechanical press was used force the 
inner hoop back into alignment with the outer hoop. Once in position, the set screws were tightened to 
establish preload on the stacks. The final prototype assembly is shown in Figure 18. 



Figure 17. — Coni 1b prototype, (a) Partial assembly, (b) Separating hoops to insert link. 



Figure 18. — Conll prototype. Shown with an aluminum mounting fixture inside. 


NAS A/TM— 20 14-218337 


14 


2.3 Damping and Power Generation Testing 

Based on the analytical work presented in the next section of this report, a method and metrics were 
developed to assess the stiffness, damping, and power generation of the vibration ring prototypes. Testing 
is being done using the high speed load frame illustrated in Figure 19(a). This system is typically used for 
elastomer testing and has the capability to apply controlled vibratory loads up to 1000 Hz. For this work, 
the system was augmented to include a data acquisition system to acquire high speed and high voltage 
sensor signals. A test fixture, shown in Figure 19(b), was built to hold the vibration ring during testing 
and rotate its orientation. The fixture includes a load cell to measure the applied force and nano-precision 
displacement probes to observe deformation. The piezo stacks are wired to external shunt circuits, which 
provide access points to measure voltage. The force, displacement, and circuit voltages are acquired by 
the data acquisition system. Preliminary data was gathered in phase 1 of the project, but it was not until 
phase 2 that systematic testing methods were established. Therefore, the performance data for all 
prototypes will be presented together in the phase 2 report. 




Displacement 

probes 


C-2013-2417 


Platen 


Q..-M 


Wired to 


external circuits 


Load cell 


Figure 19. — Test setup, (a) Load frame, (b) Test fixture. 
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3.0 Analysis 

In this section, a 2D analytical model is created to represent the dynamic electro-mechanical response 
of the vibration ring in closed form. The compression cage, piezo stack, and electric circuit features are 
related to the frequency-dependent stiffness and loss factor of the vibration ring. FEA is used to evaluate 
the analytical model assumptions and to provide greater insight about the mechanism. The analytical 
model is then used to study how the vibration ring parameters relate to damping performance using the 
previously introduced driveline model. 


3.1 Shunted Piezo Stack Model 


As a first step, the shunted piezo stack is formulated as a complex spring, so it may be incorporated 
with the mechanical framework of the vibration ring. Relationships are also derived that allow for the 
stack damping to be maximized and for power generation to be predicted. Such equations may be found 
in various forms in the technical literature (e.g., (Hagood & von Flotow, 1991)). Nonetheless a complete 
and systematic development is presented here to provide insight. For reference, the scope and 
assumptions of the model are listed in Table 2. 

TABLE 2. — SCOPE AND ASSUMPTION OF THE SHUNTED PIEZO STACK MODEL 

1 Each element of the piezo stack is modeled using the linear constitutive equations of piezoelectricity (ANSI, 1987). 

Therefore nonlinear effects are not considered. 

2 Only the motion along the stack’s poling axis (the 3-axis) is described and so numerical subscripts designating the 

input and output variable orientations have been omitted. 

3 Unless otherwise stated, waveforms are represented by phasors with the following convention. Results are therefore 
indicative of steady state response to sinusoidal type excitation. 

y (0 = y 0 cos ^ ) = Re [ye J(at } => y = y 0 eJi>y ( 5 ) 

4 Electric charge leakage and mechanical hysteretic losses are neglected. 

5 When relating properties of the piezo elements to stack properties, the following assumptions are made: 

a) Stress and charge are assumed to be evenly distributed across the face area of each element. 

b) Strain and electric field are assumed to be uniformly distributed through the thickness of each element. 

c) The compliance of the insulating end-plates and electrodes/glue between elements are neglected. 

d) The in-plane constraints created by bonding the elements are not taken into account. 

However, when the model is parameterized based on measured piezo stack data these assumptions are not required. 

Considering a single piezo element, the linear constitutive matrix equation for electro-mechanical coupling (ANSI, 
1987) is, 



( 6 ) 


The equation represents mechanical strain, S , and electric displacement, D , as functions of applied 
mechanical stress, T 7 , and electric field, E. The relationships are governed by the element’s elastic 
compliance constant, s E , permittivity constant, and its piezo constant, d. The compliance constant has 
the superscript E to designate that it was measured under a condition of constant electric field (i.e., when 
the electrodes were shorted). The permittivity constant has the superscript T to designate that it was 
measured under a condition of constant stress (i.e., with free mechanical boundaries). The piezo constant 
represents the special property of piezo materials, which allows for stress to cause electric displacement 
and electric field to cause strain. By assuming that the stress and charge are evenly distributed across the 
element’s face area and that the strain and electric field are uniformly distributed through the element 
thickness, the following change of variables is made, 


' yjl 1 = L £ d^Fa/A 

c j ,j _v a /i _ 
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Figure 20. — Sign conventions for piezo element variables. 


The output variables have become displacement, y, and current, /, while the input variables are force, F, 
and voltage, V. The geometric factors, A and /, are the element’s face area and thickness. The subscript el 
is used to distinguish the elemental quantities represented here from those associated with a complete 
stack. To clarify the sign conventions, Figure 20 shows the input and output variables with respect to an 
element. 

For a typical piezo stack, several elements are connected in a mechanical series, but with their 
electrodes connected in parallel. The combination of ^-elements into a stack is mathematically described 
by the following equations: 


F = F eb y = n y e b V = V eU 1 = nI e\ ■ 

The variables listed without subscripts are associated with the complete stack, rather than a single 
element. Using the above transformations, Equation (7) is rewritten as, 


7 


i 

S 

^ sr 

nd 

~F~ 

J/joy 


nd 

nAs T 

V 



i 



( 8 ) 


( 9 ) 


The following approximations are then introduced to relate the element properties to the stack properties, 


^sc 


A nAz 

, C = — ; — , d Qff =n-d 


ills 


l 


( 10 ) 


Here k sc is the short-circuit stiffness or the stiffness when the electrodes of the stack are shorted together 
and C is the stack capacitance when the stack ends are mechanically unconstrained. The factor d Q ff is the 
effective piezo constant of the stack, which relates voltage to displacement when the stack ends are 
mechanically unconstrained as well as relating force to charge (7//oo) when the electrodes are shorted 
together. It should be noted that the short-circuit stiffness equation does not take into account the 
compliance of the insulating end-plates or electrodes/glue between the stack elements. The bonding layers 
between elements also impose in-plane constraints that have not been considered. Therefore these 
equations must be viewed as approximations. Applying the transformations leads to the following set of 
equations. 



d&V 


( 11 ) 
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( 12 ) 


— = d eff F + CV 

J © 

These are the most basic linear equations of a piezo stack, which indicate how force or voltage gives rise 
to displacement and charge. The following equation is introduced to represent the connection of a shunt 
circuit to the electrodes. 


h 



(13) 


Here z L is the load impedance and I L is the current flowing into the load. Equation (13) is inserted into 
Equation (12) to produce expressions for generated voltage and power delivered to the load in response to 
applied force. 


y _ ^eff F 


C + \/jm L 

Pl = V 2 / z L 


(14) 

(15) 


Inserting Equation (14) into (11), and solving for Fly , leads the following equation representing the 
complex valued stiffness of the stack. 


1 dj f f 

*sc C + \/j(az L 


(16) 


The complex valued stiffness carries information about both stiffness and damping, as will be explained 
later in this formulation. As shown by the equation, its value depends on the electrical impedance. The 
minimum value is the short circuit stiffness £ sc , which corresponds to load impedance z L = 0. The 
maximum value, called the open circuit stiffness, corresponds to z L = oo, and is defined by, 


^oc 


^eff 

c 


(17) 


The electromechanical coupling coefficient, k, is now introduced to simplify the complex stiffness 
equation. This variable was defined in (ANSI, 1987) as a dimensionless material property, relating the 
electrical energy generated to the mechanical energy applied. This definition is extended to a piezo stack 
by examining the conceptual energy conversion experiment illustrated by Figure 21. The experiment has 
three phases: (1) Force is applied in absence of an electrical load, requiring the applied energy to be 
stored within the stack. (2) The displacement is fixed, while a low impedance electrical load is connected 
and the electrical energy is drained. Due to the piezo effect, this causes a reduction in force magnitude. 
(3) With the electrical load in place, the piezo is returned to zero displacement. 
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Assuming no losses within the piezo stack, the electrical energy delivered to the load is W\ and the 
applied mechanical energy is W\ + W 2 . The coupling factor, defined in terms of energy, is therefore 
related to the open and short circuit stiffness quantities as shown below. 


K 


2 .. 


w { + w 2 


oc-^sc) 


2 

Tmax 



(18) 


Based on this equation, measurements of short and open circuit stiffness may be used to indirectly find 
the coupling coefficient. The more common expression for the coupling coefficient, shown below, is 
found by inserting Equation (17) into (18). 


d. 


K 


effj_sc_ 

c 


(19) 


Inserting Equation (19) into (16) leads to the following expression for complex stack stiffness. 

2 


k = k sc H, H = 1 + 


K 


2 1/ ycoC 

1 — K Z H — — 


( 20 ) 


Note that Equation (20) was previously derived and presented in (Davis & Lesieutre, 2000). With this 
representation it’s clear that complex stiffness is controlled by the impedance ratio between the piezo 
capacitance and the load. To separate the complex valued stiffness into lossless stiffness, k , and loss 
factor, r|, the following expansion is applied. 


k = k[l + jr}], (21 

The derivation of this relationship, originally presented by (Ungar & Kerwin Jr., 1962), is examined here 
to clarify the meaning of k and r|. Consider Hook’s law for a complex valued spring, 


F = ky = k[l + jr\]y. 


( 22 ) 
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The following is the equivalent expression in the time domain, 


F(t) 


,y(0 _1 D ; o sin(co? + (|) 3; ). 


( 23 ) 


By applying Pythagorean’s trigonometric identify this becomes, 

F(t) 


■=y(t) ±r \^-y^ cos 2 (©i + by )• 


(24) 


Inserting the phasor definition of y(t) from Equation (5) leads to the following relationship between force 
and displacement. 




(25) 


The equation shows that when r| = 0 the relationship between force and displacement is that of an 
ordinary spring with spring constant k. Thus, k is the lossless component of the complex stiffness. Plots of 
Equation (25) in Figure 22 illustrate that finite values of r| create a hysteresis loop in the relationship 
between force and displacement. Considering steady-state oscillation, the energy dissipated per cycle may 
be calculated based on the energy input per cycle as follows. 

w d =$F\y{t)]dy 

y 0 y Q 

= k j y{t)+r\yjyo - y 2 ( t)dy -k J y(t) -r\ ^y 2 0 - y 2 (t)dy @6) 

-yo -y a 

= Kr\kyy 


The total oscillation energy is calculated as the energy stored in the lossless component of the complex 
spring at maximum displacement (at that point kinetic energy is zero), 


W«=\Fo- 


(27) 


By combining Equations (26) and (27) the loss factor is defined as, 


fi = 


^tot 


(28) 


Thus the loss factor is a damping parameter, which is proportional to the fraction of energy dissipated 
within the complex spring during a cycle of steady-state oscillation. 

Applying the expansion (21) to the complex stiffness expression (20), the loss factor for the piezo stack 
is found to be, 


*1 = 


9 

K CO Cz L 

l + (l-K 2 )(co Cz L f 


(29) 
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Figure 22. — Example force vs. displacement cycles based on Equation (25). 


The expression shows that the loss factor is frequency dependent and its magnitude is affected by choice 
of load impedance. For the limited case of resistive loads, the equation dr\/dR L = 0 is solved to find the 
value that maximizes the loss factor, 


D max 

K L 


1 

COCVl-K 2 


(30) 


By inserting Equation (30) back into Equation (29), the maximum value for the loss factor is expressed as, 


0 


max 



(31) 


This equation shows that the maximum loss factor is completely determined by the coupling coefficient. 

In the design work implemented using this stack model, Equations (30) and (31) are used to select the 
stack and associated load resistance values that maximize the loss factor in specific frequency bands. 

Then Equation (20) is used to represent the stack as a complex spring element within the vibration ring 
mechanical framework. The system is simulated to produce the steady state vibration response of the 
vibration ring, and the corresponding forces transferred to the stack. Finally, the stack forces are used to 
calculate the current, voltage, and dissipated power using Equations (13) to (15). 

There is a set of five stack parameters used in the model, {ic,<7 eff ,C,£ sc ,£ oc } . Three of these should 

be determined from data, and the remaining two computed based on the linear relationships in 
Equations (18) and (19). As indicated by Table 3, the parameters may be determined based on element 
properties or stack measurements. Element properties are often available from manufacturer data sheets 
making them convenient to use. However, the equations relating element properties to stack parameters 
require several simplifying assumptions (those listed in Table 2, row 5). It’s therefore preferable to use 
direct stack measurements to derive the model parameters. Table 3 shows the required stack 
measurements, though standard data collection methods do not exist at this time. 
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TABLE 3.— OPTIONS FOR DETERMINING THE STACK MODEL PARAMETERS 



Stack model parameters | 

Via 

K 

<Cff 

C 

k sc 

k oc 

Element properties 21 

K el 

n-d 

CO 

A 

nls E 


Stack 

measurements 21 


Ay A 7/ z co 

AF F=0 A F v=0 

A//yco| 

AF 

A f\ 

W x +W 2 

> 

II 

o 

A y z L =o 

K 

II 

8 


Parameters based on element properties require the extra set of assumptions listed in Table 2, row 5. It’s therefore more accurate 


to determine parameters from stack measurements. 


TABLE 4.— SCOPE AND ASSUMPTIONS OF THE VIBRATION RING MODEL 


1 

Waveforms are represented by phasors to be consistent with the stack model. 

2 

Only radial force excitation is considered. 

3 

Deformation is assumed to be local, such that a 2D model represents the response of a finite circumferential section of the 
vibration ring. 

4 

Kinematic motion of the compression cage is modeled using rigid links (fixed length) that rotate about pin joints. 

5 

The effective cage stiffness, 7r C a g e, is assumed to be linear. Linearity permits the use of superposition to incorporate the 
stack stiffness, k , into the compression cage framework. 

6 

The effective cage stiffness in the y-direction is assumed to be negligible in comparison to the piezo stack stiffness, k . 

7 

The compression cage stiffness ratio, a, is assumed to be independent of the link angle, so that it can be used as an input 
variable. 

8 

Inertial forces are not considered. 


3.2 Vibration Ring Model 

Using the complex spring representation of a shunted piezo stack, developed in the previous section, a 
2D mechanical model of the vibration ring is developed. First, kinematic relationships of the compression 
cage are written, relating external compression to compression of the piezo stack. Springs are then 
introduced to represent the bending resistance of the compression cage as well as the complex stack 
stiffness. Free body analysis leads to a system of equations to predict the response of the vibration ring to 
vibratory excitation. These equations are then linearized to represent the vibration ring by its stiffness and 
loss factor parameters. These parameters are defined in relationship to the stiffness and loss factor of the 
piezo stack. The model is developed using the scope and assumptions listed in Table 4. 

Figure 23 shows a cross sectional drawing of one of the more mature vibration ring design concepts 
(Coni la), both in its unloaded and loaded state. In this design, the compression cage is made of several 
butterfly amplifier segments, each paired with a discrete piezo stack. It should be noted that the model 
developed here may be applied to a design that uses a single ring-shaped stack by assuming that each 

amplifier segment is supported by a section of the stack circumference with associated stiffness, k . The 
figure has line segments superimposed to show the important dimensions of the compression cage. The 
geometric variables that change under load are shown with subscripts ‘unload’ and Toad ’ to indicate their 
state. The link and wall lengths do no change under load, and are defined by the following two equations. 


t _ ^unload ^ 
Mink 0 . 

2cosc|) unload 


(32) 


Avail Tunload + (-^unload w ) ^ an ^unload 


(33) 
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(a) Unloaded geometry 


(b) Loaded geometry 



Figure 23. — Vibration ring model geometry based on Coni la. (a) Unloaded, (b) Loaded. 


The motions of the compression cage are described by the following set of equations. The equations 
indicate that horizontal compression, x, increases the link angle, cp, which in-turn increases vertical 
compression, y'. 1 


^load -^unload x 


(34) 



f \ 



4>load = 

= cos 1 

*load- W 
V ^Aink J 


(35) 

= Avail 

~h L 

link ( "^load 

- w'f 

(36) 

y 

y unload T load 


(37) 


Figure 24 shows a simplified dynamic model (Fig. 24(a)) and corresponding free body diagrams 

(Fig. 24(b) and (c)). The complex stiffness, k , is included to represent either a shunted piezo stack or the 
load section of a shunted ring-shaped stack. In the simplified model, the walls and links are rigid and 
interconnected using pin joints. In reality the joints resist bending, which creates an effective stiffness in 
the horizontal and vertical directions. Here the vertical stiffness is neglected, because it is assumed to be 
small compared to the stack stiffness. The effective horizontal stiffness is captured by adding linear 
springs at the top and bottom, with combined stiffness £ cage . This value is also meant to include any other 
sources of cage resilience, such as the bending stiffness of the outer hoop. Based on the free body 
diagrams, the forces that arise due to compression of the structure are defined as follows. 


i 


The prime reverses the sign convention. Primed variables therefore associate a positive sign with compression. 
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(a) Model 


(b) Free body- outer wall (c) Free body- top plate 



Figure 24. — Dynamic representation of vibration ring, (a) Lumped mechanical model, (b) Free body 
diagram of outer wall, (c) Free body diagram of stack top plate. 


F =k-y 


F = k . r 
± cage 'xage 


^link “ ' 


sin (|>i oac i 

Summing the forces on the outer wall provides the following equation for the external force. 

^x,=^cage+A* C 0 S <L d 


(38) 

(39) 

(40) 


(41) 


The set of equations will now be combined and linearized to produce expressions for the complex 
vibration ring stiffness as well as the transfer of external force to the stack. Combining Equations (32) to 
(41) provides the following expression for the complex stiffness, 


p 

^vr — ^cage ^ ^ (^unload ) 


(42) 


The function G, expanded below, relates the complex stack stiffness to the complex vibration ring 
stiffness. 


c (^unload ~X~w) 1 
X 


( -^unload ^ 

/ (^unload _ 

V COs2 ^unload 


) ^unload 


(•^unload 



(43) 


Strictly speaking the value of G is dependant on the deformation response, x, and so it’s a nonlinear term. 
However, by assuming that the value of G changes negligably with x, the following linear approximation 
may be used. Note that this approximation was derrived by evaluating G at x = 0, applying L’hopital’s 
rule, and limiting the domain to (0,90) degrees. 
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Q _ | ! CQS ^ ( l ) unload 
Sin2 Unload 


( 44 ) 


Figure 25(a) shows the exact G function plotted over a practical range of deformation and for several 
different unloaded link angles. The function at x = 0 represents the estimate. The plot shows that the value 
of G deviates only slightly as x is increased. Figure 25(b) shows the complex stiffness using the estimate 
of G. Note that the results are shown normalized to the exact value of complex stiffness, and that £ cage has 
been set to zero as a worst case condition. The plots shows that stiffness is overestimated by less than 
1 percent, except when the link angle is lower than 35°. 

Using the linear approximation for G, the vibration ring is treated as a complex- valued linear spring 
and expanded into its lossless stiffness component and loss factor according to Equation (2). Using this 
expansion as well as the equivelant stack stiffness expansion (21), the expression for complex vibration 
ring stiffness (42) is rewritten as, 

^vr ( 1 "I" /Hyr ) — ^cage + k ( 1 + /T|) • G ( (|) un i oa d ) • (45) 

By rearranging the real part of this equation, the following expression for stiffness is found. 


k VY k'G( ^ un i oac [ ) + ^cage 


Then by defining the stiffness ratio, 


(46) 


« = W 

Kr 

the equation for vibration ring stiffness is simplified to, 

, k ’ G (^unload) 


(47) 


(48) 


By rearranging the imaginary part of Equation (45) the expression for the vibration ring loss factor is 
found to be, 


fivr 


^unload) 1 ! 


(49) 


Inserting Equation (48) reduces the expression to, 

rivr = (l — a ) T l- ( 5 , 

Equations (48) and (50) show the basic relationships between the physical properties of the shunted 
stack and vibration ring. The vibration ring stiffness is proportional to the stack stiffness, with the 
proportionality constant (G) controlled by the link angle. As illustrated in Figure 26(a), different 
combinations of stack stiffness and link angle can be selected to produce the same target vibration ring 
stiffness. At low link angles, there is a stiffness magnification effect and therefore a very soft stack may 
be used. At high link angles, the opposite is true. Moreover, the vibration ring loss factor is proportional 
to that of the stack, and is not affected by the choice of link angle. Therefore, as illustrated in 
Figure 26(b), each of the stack-angle combinations that produce the same vibration ring stiffness has the 
same FT. The compression cage stiffness, as defined by a, has the effect of increasing the vibration ring 
stiffness, while decreasing the loss factor. It should be minimized to achieve maximize damping 
performance. 
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Figure 25. — Plots supporting linear approximation of G. (a) The G-parameter is shown to vary only slightly. 

(b) The resulting vibration ring stiffness is overestimated by less than 1 percent, except when the link angle 
is lower than 35°. 



^ unload ( de 9) 



Frequency (Hz) 


Figure 26. — Adjustability of the vibration ring, (a) Plot of link angle and stack stiffness combinations that produce 
the same vibration ring stiffness, (b) Simulation of simplified driveline (Eqs. (1) to (4)) showing that each 
combination produces the same force transmissibility. 


2 The relationship shown is for constant stack stiffness, but the same idea applies when the stack stiffness 
frequency dependent. 
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The following expression for compressive stack force is found by combining Equations (39) to (42), 


F 


Fqx t ^ an ^load 


f 

1 - 

V 


k \ 

^cage 

b 

V vr J 


(51) 


This equation is also non-linear, as it depends on the response angle cpi oa d- The loaded link angle is 
replaced by the unloaded link angle to arrive at the following deterministic approximation. 


F ^ext ^unload 


f k \ 

j 'xage 

K, 


V vr y 


(52) 


Figure 27(a) shows that the exact and estimated stack forces are very similar. Here k cag e = 0, as this 
produced the greatest deviation. Figure 27(b) shows the estimated stack force normalized to the exact 
value to better quantify the deviation. The data shows that stack force is underestimated by about 
1 percent or less. 

Rearranging the stack force estimate equation produces the force transfer ratio metric, 


FTR = 


^ext ^unload 


= 1 


v cage 


(53) 


This equation is used to compare efficiency of transferring force to the stack. The equation also reveals 
that the cage stiffness essentially acts as a force shunt around the stack, which explains why reducing a 
increases the vibration ring loss factor. 




Figure 27. — Plots supporting linear approximation of F. (a) The exact and estimated values are shown to be similar 
(b) Stack force is underestimated by about 1 percent or less. 
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Using the linearized equations, the vibration ring performance is completely defined by the stack 
parameters and the compression cage variables a and cp un ioad. None of the other geometric variables are 
needed. The equations are used in two different linear combinations. In the next section, the stiffness and 
FTR equations are used to compare with static FEA results for ConlO. Subsequently, the stack and model 
equations are used dynamically to determine how specific parameters influence damping. 


3.3 Vibration Ring Finite Element Analysis 

As a part of the design process, FEA models were developed for ConlO and 11. This was done to tune 
the stress distributions and optimize the FTR metric (53). As an example, Figure 28 shows the stress 
distribution from Coni Ob. In this section, the FEA model is examined in more detail and simulation 
results are compared with analytical model predictions. 

The features and limitations of the analysis are listed in Table 5. 



0.000 


0.500 



2.000 (in) 


Figure 28. — Example FEA results showing stress distribution for Con 10b. 


TABLE 5. — FEATURES AND LIMITATIONS OF THE FINITE ELEMENT ANALYSIS 


1 

Only static analysis was implemented. 

2 

Stacks are represented by springs, and the spring force is spread over the contact area of the stack face. 

3 

Stack stiffness is real-valued, because loss factor is always zero for static loads. 

4 

Compression cage material properties are for stainless steel. 

5 

Assembly screws in ConlO are not represented explicitly. The mating surfaces are constrained to have 
matching displacements, moving as if bonded together. 

6 

Force is applied to the outer hoop with a flat rigid surface, as illustrated in Figure 28. 

7 

The inner hoop has fixed nodes so that it remains round. 

8 

Specific information regarding the prototypes is provided in the design section of this report. 
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Parametric study was conducted using the ConlOb model following the steps in Table 6. Figure 29 
provides further insight into how the model was tuned at each link angle. Given that the hoop diameters for Con 10 
and 1 1 were limited by the design specifications, changes in link angle were accommodated by changes in the link 
length. I.e. links were made to be longer when the angle was increased. As indicated by the procedure in Table 6, 
each design was tuned under a 400 N load by adjusting the radii and thickness dimensions shown in Figure 29 (Rl-4 
and Tl-2). The first objective was to create even stack compression. To accomplish this, the inner links needed to be 
thinner and have smaller radii than the outer links. The second objective was to maximize the average stack 
compression, which is equivalent to maximizing the FTR as defined by Equation (53). This was a balance between 
establishing a rigid path for the force to transfer longitudinally through the links (resulting in high k vr ), while 
allowing the links to bend with little resistance (resulting in low k cage ). Clearly these were opposing goals with the 
present design concepts. The final objective was to minimize the maximum stress. Maximum stress was generally 
located within the link radii where the most bending occurred. 

Static analysis was done with the analytical model for comparison with the FEA results. The 
calculation steps are listed in Table 7. Here the k CSL ge values were taken from the FEA results in order to 
highlight the other differences that exist between the models. 


TABLE 6.— FEA PROCEDURE FOR PARAMETRIC STUDY OF CONI OB 


1 

Select Unload 

2 

Set F ext = 400N 

3 

Adjust the link geometry (radii Rl-4 & thicknesses Tl-2 from Figure 29) to achieve: 

• Even stack compression 

• Maximum average stack compression 

• Lowest maximum stress 

4 

Apply static forces up to 845 N and document x,F , and maximum stress.* 
*Note that stack force was only recorded for 45° 

5 

Remove the stack springs and repeat step 4. 

6 

Compute the following 

W =— > =— . FTR= — — 

x stacks removed X Uxt t unload 


compression 

Figure 29. — Graphic to explain link geometry tuning procedure. 




TABLE 7.— ANALYTICAL MODEL CALCULATION STEPS FOR COMPARISON WITH FEA 


1 

Input static parameters 

l W F; xt ,(t) unload 

constant fr^raA Iterate 

2 

Compute vibration ring stiffness 

G = 1 + co S 2V l0ad tK=k ^ +kG 

Sm vunload 

3 

Compute stack force transfer ratio 

FTR = 1 kca °° e 

^cage 
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Normalized max vibration ring stress vs. deflection 



Deflection (m) 


-♦-40 

-*-50 

-*-55 


Figure 30. — FEA prediction of normalized maximum vibration ring stress vs. deflection. 
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3.3.1 Stress Analysis 

The maximum stress associated with each design angle and loading condition is plotted in Figure 30 
against deflection. The stress is normalized to deflection to allow for comparison between designs that 
have different stiffness. The flatness of the curves indicates that stress is a linear function of deflection for 
each design angle. Also, higher angle designs have lower normalized stress. This can be attributed to the 
increase in link length. For the same deflection, a longer link has less strain and therefore lower stress. 

3.3.2 Evaluation of Analytical Model Assumptions 

Several assumptions were made to develop the 2D analytical model. Each of the assumptions that 
relate to the mechanical response is evaluated here by examining the FEA results. 

Assumption #3: Deformation is local 

The assumption of local deformation was made so that the 2D model only needed to represent a 
single circumferential section of the vibration ring. For this to be valid, the 3D vibration ring response 
must only deform and experience stress in the vicinity of the applied load. Figure 31 shows exaggerated 
deformation patterns and stress distributions for ConlOb and Coni lb. Referring to Figure 31(a), most of 
the deformation in Con 10 appears to be local. However, there is some noticeable hoop deformation 
outside of the loaded region. In addition, there is fairly significant stress in the links next to the loaded 
section. The wide mounting tabs appear to be increasing the bending stiffness of the hoop, causing the 
applied load to be distributed circumferentially. As shown in Figure 31(b), the results for Coni 1 are 
better. There is relatively little hoop deformation outside of the loaded region and minimal stress in the 
neighboring links. Despite the differences between these designs, the assumption of local deformation is 
considered to be reasonable. 
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Figure 31. — Deformation patterns and stress distributions in vibration ring FEA models, (a) ConlOb, (b) Conllb. 



Figure 32. — Deformation patterns and stress distribution in links, (a) ConlOb, (b) Conllb. 


Assumption # 4 : Links are rigid 

In the analytical model, the links were assumed to bend at the hinge locations but remain rigid in- 
between. Plots of the FEA deformation and stress patterns are shown in Figure 32. The link pattern of 
ConlOb, shown in Figure 32(a), is shaped like an ‘S’. Rather than bending only at the hinge locations it 
bends globally. There is less deformation within the links of Coni lb, as shown in Figure 32(b), but the 
links pivot at the top. Neither deformation pattern completely matches the assumption. 
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Given this discrepancy, there was a concern that the links may be softening the force transfer path 
between the vibration ring and stacks. FEA simulation was done with Con 10b and Coni lb where the 
stacks were made to be completely rigid. This enabled a measurement of stiffness, called blocking 
stiffness, which was associated with the link deformation. As shown in Table 8, the blocking stiffness was 
found to be on the same order of magnitude as the stack stiffness. This suggests that the links have a 
significant softening effect on the mechanism, which is not represented by the analytical model. 

Assumption #5: Cage stiffness is linear 

The linearity of £ cag e is assumed when the stack spring, k , is inserted into the compression cage 
framework. This is an implicit use of the principal of superposition for linear systems. As a check of this 
assumption, the FEA model for Con 10b was exercised with the stacks removed and the force and 
deflection were documented. The resulting cage stiffness curves for each design angle are shown in 
Figure 33 as a function of the applied force. It’s apparent from the data that cage stiffness is quite linear. 

Assumption # 6 : Alpha is independent of link angle 

The parameter a was assumed to be independent of link angle, so it could be used as an input 
parameter in the analytical model. Recall that a is the ratio of cage stiffness to vibration ring stiffness. To 
evaluate the assumption, the FEA cage stiffness, vibration ring stiffness, and their ratio are shown in 
Figure 34 verses link angle for ConlOb. Several loading conditions from 44 to 845 N are plotted here, but 
for the most part the curves lay on top of each other. Cage stiffness and vibration ring stiffness both 
decrease with link angle. However, they decrease at similar rates and their ratio remains constant at 
approximately 0.3. This suggests that the assumption is valid. 


TABLE 8.— BLOCKING STIFFNESS 
COMPARED TO STACK STIFFNESS 
(Data collected from 45° models at 400 N) 



Stack 

(N/jam) 

Blocking 

(N/jLim) 

ConlOb 

68 

49 

Coni lb 

60 
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6.00E+06 
4.00E+06 
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0.00E+00 h 1 1 1 1 1 1 1 1 1 

0 100 200 300 400 500 600 700 800 900 

Applied force (N) 

Figure 33. — FEA prediction of cage stiffness vs. force. 
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Cage / Vibration ring stiffness Stiffness (N/m) Stiffness (N/m) 
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1.40E+07 

1.20E+07 

1.00E+07 
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6.00E+06 

4.00E+06 

2.00E+06 

0.00E+00 


(a) Cage stiffness vs. link angle 



t 1 1 1 1 1 r 
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(b) Vibration ring stiffness vs. link angle 
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(c) Stiffness ratio vs. link angle 
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0.00 -I T T T T T T T 
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Figure 34. — FEA prediction of stiffness vs. link angle, (a) Cage stiffness, (b) Vibration ring stiffness, 
(c) Stiffness ratio. 
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3.3.3 Comparison With Analytical Model 

Here the FEA predictions of vibration ring stiffness and FTR for the Coni Ob design are compared 
with the analytical model. Vibration ring stiffness versus link angle is plotted for all load cases in 
Figure 35. The bottom set of curves are the FEA results and the top set are the analytical results. The data 
shows that the analytical model significantly over-predicts the vibration ring stiffness, especially at the 
lower angles. Force transfer ratio data was only collected at 45°, and so it is plotted against applied force 
in Figure 36. The data shows that the analytical model also overestimates the FTR. It is hypothesized that 
these discrepancies are due to the assumption of rigid links being incorrectly applied in the analytical 
model. There analytical model should therefore be modified to represent the finite link stiffness. 



Figure 35. — FEA and analytical model predictions of vibration ring stiffness. 


Force transfer ratio at 45° 



Applied force (N) 

Figure 36. — FEA and analytical model predictions of force transfer ratio. 
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3.4 Driveline Damping Study Using Analytical Model 

With an understanding of the strengths and limitations, the analytical model is used to study how the 
Vibration ring configuration affects driveline damping. The equations indicate that the parameter k 
controls the peak loss factor of the stack (Eq. (31)), while a determines how much of the stack loss factor 
is reflected to the vibration ring (Eq. (50)). To understand the relative influence of these damping related 
parameters, a numerical study is done using the driveline model of Figure 2 (Eqs. (1) to (4)). Several 
vibration ring configurations are created by iterating k and a though realistic ranges. Each configuration 
is simulated in the driveline model to compare damping performance. The calculation steps are listed in 
Table 10, and summarized in the text below. 

The process starts by assigning several input parameters, including a, k, the link angle, and p. The 
link angle is set to 45°, but as discussed in the previous section, from a conceptual design perspective this 
has no influence on the properties or performance of the vibration ring. Of course, selection of the link 
angle will influence details of the design. For example, it changes the stack stiffness needed to achieve a 
specific value for vibration ring stiffness. 


TABLE 10.— CALCULATION STEPS FOR DRIVELINE DAMPING STUDY" 


1 

Input 

^unload /^^OC 
Constant Iterate 

2 

Calculate system information 

m s = . b ,2> k b = k b( l + J r \b) 

(2 *f nJb ) 

3 

Set vibration ring stiffness to 
create specific frequency shift. 

fn = P. fn,b > K«,n = ™ s ( ) , K,n = , ’ 

K b ^eff ,n 

4 

Set RC combination to 
maximize loss factor at new 
natural frequency. 

1 K 2 

R'TC- , T| max - * ,<. ax -(l cx) r| max 

2nf n ^\ - k 2 2V1-K 2 

5 

Relate vibration ring stiffness 
to stack stiffness at new natural 
frequency. 

Q_ \ \ C0S ^unload £ _ _ a )Kv,n 

Sin 2 <t>unload ’ ” G 

Hn ” 1 + 1 - k 2 + Vj2nf„RrC ’ ^ ~ Re {H n } 

6 

Compute complex stiffness 
spectra 

r K 2 1 ~ 

cage vr - n ’ L l-K 2 +1/ j2nfR^ x c\’ k sc 

K = *cage + kG, Cff = - b 7 

h+K T 

7 

Compute driveline 
performance spectra 

Dynamic compliance: c = — ^ ^ — 

^in ~(2j if) m s +k eff 

Force transmissibility: FT = c • k eff 

8 

Compute forced response 
spectra 

Force on vibration ring: F ext = F out = F m • FT 
Stack force: F = F ext tan ()) unload [ l - 8 

1 K ) 

k 2 (-F') 2 

Electrical power: P L = — 

**(1+1/ nmipc) Rrc 


"Variables defined at the natural frequency have the subscript n 
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The P parameter defines the natural frequency reduction that the driveline system can tolerate when 
the vibration ring is integrated. Reduction in natural frequency is associated with softening of the 
vibration ring. This improves the system damping performance, but the allowable shift is generally 
constrained by the operating frequencies of the machine. Therefore the selection of P in practice will be 
system specific. Here a value of 0.8 (20 percent reduction) is arbitrarily selected for all simulation cases. 

Using the value for p, the vibration ring stiffness is chosen using the formula for un-damped natural 
frequency in a single degree of freedom system. Then, the RC combination of the circuit and stack is 
selected with the purpose of maximizing damping at the new natural frequency. Next, the required short- 
circuit stiffness of the stack is found. Finally, the complete spectra of stack and vibration ring properties 
as well as driveline performance metrics are calculated. 

3.4.4 Case 1: Variation of k (piezo electromechanical coupling factor) 

For simulation case 1, k is iterated through the set {0.5, 0.7, 0.9} while a is fixed. Note that 0.5 
represents a very weak electromechanical coupling factor, 0.7 is typical of piezoceramic stacks, and 0.9 is 
typical for single-crystal piezo stacks. Figure 37 shows the spectra of stack and vibration ring properties. 
Here the dashed vertical line indicates the predicted new natural frequency after the vibration ring is 
integrated with the driveline. As intended, the vibration ring stiffness for each case is identical at the 
natural frequency. However, the frequency dependence of vibration ring stiffness increases with k. For 
example, the configuration with largest k has both the lowest DC stiffness and the highest stiffness at the 
maximum frequency. The stack stiffness follows the same trends as the vibration ring stiffness, but its 
magnitude is slightly smaller due to the finite value of a. 


250 

200 

IT 

n. 

r 150 
100 





f/fn(-) 


Wn(-) 


Figure 37. — Influence of xon property spectra (simulation case 1). (a) Vibration ring 
stiffness, (b) Stack stiffness, (c) Stack loss factor, (d) Vibration ring loss factor. 
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As expected, the stack loss factor peaks at the natural frequency and drops to zero at DC. The 
vibration ring loss factor follows the same trends, though the absolute values are 80 percent less due to 
the finite value of a. Changes in k appear to have a profound effect on the vibration ring loss factor. For 
example the change from k = 0.7 to 0.9 causes a 270 percent increase in the loss factor peak. This drastic 
improvement motivates the use of single-crystal piezo stacks over piezoceramic stacks. 

The driveline system performance for simulation 1 is shown in Figure 38. In the FT plot, the natural 
frequency was shifted to the target location for the lower two cases of k. However, at the highest value of 
k the natural frequency is slightly greater. This is because the calculations were based on the formula for 
un-damped natural frequency, which is imprecise for heavily damped systems. The calculation procedure 
may be considered conservative, since the natural frequency was shifted to a larger value than planned. 


Force transmissibility 




Figure 38. — Influence of k on driveline system performance spectra 
(simulation case 1). (a) Force transmissibility, (b) Dynamic compliance, 
(c) Electrical power. 
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For all cases the reduction in peak FT response was significant, on the order of 10’s of dB. Note that the 
exact reduction in peak FT is not relevant, because the driveline damping parameters (with no vibration 
ring) were selected arbitrarily. The variation in peak FT between the cases was significant, 13 dB, 
supporting the notion that k is a critical parameter. Reduction in peak FT is directly related to the peak 
loss factor, indicating that loss factor is an appropriate metric for evaluating the vibration ring’s potential 
damping performance. 

The trends in dynamic compliance around the natural frequency are nearly identical to the FT, which 
indicates that a vibration ring design that reduces transmitted force will concurrently suppress the motion 
of the shaft at resonance. On the other hand, the trends for power generation are unique. Below and above 
the resonance the power generated is directly related to k. This is expected, because stacks with higher k 
convert more of the applied mechanical energy to electrical energy. Around the resonant peak, however, 
the trend reverses; as k is increased less power is generated. I.e. when the resonance is dampened, less 
work is done on the vibration ring and therefore less power is generated. To make use of the generated 
energy, of course, the resistive load would need to be replaced by an energy harvesting circuit with the 
desired impedance. 

3.4.5 Case 2: Variation in a (Stiffness Ratio) 

For simulation case 2, the influence of cage stiffness is studied. The ratio of cage stiffness to overall 
vibration ring stiffness, parameter a, is iterated through the set {0.3,0.2,0.1} while k is fixed. The range 
for a was selected based on the intuition gained after creating several compression cage designs. The 
highest value is meant to correspond to a preliminary design and the lowest value represents an optimized 
design. Note that for a compression cage that transmits significant torque, the value of a may be higher 
than 0.3. Figure 39 shows the spectra of stack and vibration ring properties, while Figure 40 shows the 
driveline performance spectra. Note that the green line in these plots represents the same configuration 
(k = 0.7, a = 0.2) as the green line in the simulation case 1 plots. 

Considering first the vibration ring stiffness, the change in frequency dependence for the different 
cases is relatively small and therefore the variation in DC stiffness is minor. The stack stiffness curves are 
shifted vertically according to the differences in cage stiffness with each design. I.e. when a is reduced it 
requires a slightly higher stack stiffness to achieve the target vibration ring stiffness. For all three cases, 
the stack loss factor spectra are the same. This is because each case has the same k (stack type) and is 
optimized to the same frequency. However, higher values of a lead to lower vibration ring loss factor. As 
mentioned previously, the cage stiffness acts as a force shunt reducing the influence of the stack. The 
resulting changes in vibration ring loss factor are significant, but much less than the changes observed 
through variation of k. Therefore, a should be considered a secondary design factor with regard to 
damping. 

All of the trends in driveline performance, shown in Figure 40, are consistent with simulation 1. 
Increases in loss factor, creating by reducing alpha, correspond to reductions in the peak FT. Here, the 
peak FT only varies by 2 dB as the variation in loss factor is minor. The dynamic compliance around the 
resonance follows the same trends as FT, supporting the idea that the vibration ring simultaneously 
reduces transmitted force and shaft motion. Considering the power spectra, the vibration ring with the 
best damping performance generates the most power above and below resonance, but the trend reverses 
close to the peak. 
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Figure 39. — influence of a on property spectra (simulation case 2). (a) Vibration ring stiffness, 
(b) Stack stiffness, (c) Stack loss factor, (d) Vibration ring loss factor. 


NASA/TM— 20 14-2 18337 


39 


|PJ (W) 20l °9lol X s /F inl < dB > 


Force transmissability 



Electrical power spectrum 



0 0.5 1 1.5 


Figure 40. — Effect of a on driveline system performance spectra (simulation case 2). (a) Force 
transmissibility, (b) Dynamic compliance, (c) Electrical power. 
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4.0 Conclusion 

4.1 Summary 

Introduction : The vibration ring was conceived as a driveline damping device to prevent structure- 
borne noise in machines. It has the appearance of a metal ring, and can be installed between any two 
driveline components like an ordinary mechanical spacer. Damping is achieved using a ring-shaped 
piezoelectric stack that is poled in the axial direction and connected to an electrical shunt circuit. 
Surrounding the stack is a metal structure, called the compression cage, which squeezes the stack along its 
poled axis when excited by radial driveline forces. The stack in turn generates electrical energy, which is 
either dissipated or harvested using the shunt circuit. Removing energy from the system creates a net 
damping effect. The vibration ring is much stiffer than traditional damping devices, which allows it to be 
used in a driveline without disrupting normal operation. In phase 1 of this NASA Seedling Fund project, a 
combination of design and analysis was used to examine feasibility of this concept. 

Design : Several design concepts were evaluated using solid modeling, finite element analysis, and by 
creating prototype hardware. Early concepts focused on developing a compression cage structure that 
could transfer radial force through the poled axis of the stack. The so called ‘segmented butterfly 
amplifier’ was found to provide the required kinematics. It used several mechanical amplifiers spaced 
around the circumference, which squeezed the piezo stack in the axial direction when the walls of the 
compression cage were moved closer together. Complete prototypes were made, using both additive and 
traditional manufacturing techniques. Additive designs had fewer parts and more straightforward 
assembly procedures, but the parts lacked the required precision. Designs based on traditional 
manufacturing were more precise, but they required several small parts and were cumbersome to 
assemble. To evaluate the damping properties and resulting power generation of the prototypes a load- 
frame test was developed. Test results will be reported in phase 2 of the project. 

Analysis’. Analytical modeling was used to evaluate the combined electro-mechanical response of a 
simplified cross-section of the vibration ring, while finite element analysis served to provide greater 
insight into the mechanical response of the detailed three-dimensional designs. In the analytical model the 
shunted stack was formulated as a complex-valued mechanical spring integrated within the mechanical 
framework of the butterfly amplifier. The model revealed that the stiffness and damping of the stack and 
overall vibration ring were related through the amplifier’s link angle and cage stiffness ratio parameters. 
Finite element analysis showed that the detailed designs did not exhibit ideal performance, because the 
amplifier links deformed under load. This reduced the mechanical coupling between the stack and 
vibration ring, thereby causing reductions in stiffness and loss factor of the mechanism. After comparison 
with the finite element model, the analytical model was used to evaluate how the vibration ring 
configuration affected damping in a simple driveline. The electromechanical coupling factor of the stack 
was shown to have an enormous influence on loss factor and the resulting driveline damping 
performance. The cage stiffness ratio of the vibration ring was somewhat influential, but much less than 
the coupling factor. Above and below the system resonance, power generation was shown to be directly 
rated to damping. However, in the vicinity of the resonance less power was generated when damping was 
increased. 
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4.2 Lessons 

The following lessons should be taken into consideration for continued development of this 
mechanism. 

1 . Considering the simulation results, the vibration ring should be stiffer than the other driveline 
components to prevent significant shifts in natural frequency, while the vibration ring loss factor 
should be made as high as possible to reduce peak force transmissibility. 

2. The stiffness and loss factor of the vibration ring are proportional to the respective stiffness and 
loss factor of the shunted stack within. The compression cage link angle scales the stiffness, but 
does not affect the loss factor. The peak loss factor is mainly controlled by the electromechanical 
coupling factor of the stack, and the frequency of the peak is tuned by selecting the circuit 
impedance. 

3. Configuring the vibration ring to reduce force transmissibility will also suppress shaft motion 
around resonance. However, as the resonant response is suppressed less power will be generated. 

4. The segmented butterfly amplifier provides the required kinematic relationship, but has the 
following design challenges. 

a. Its geometry requires significant axial width. Therefore future designs should compensate 
for this by using shorter stacks or by modifying the amplifier geometry. 

b. To effectively transfer force through the stack, the links must be stiff, and the joints must 
bend freely. An alternative link design should be considered, which allows for the joints 
and links to be tuned separately. 

5. Assembling the compression cage around a ring-shaped stack is not trivial. Concepts 8, 10a, and 

1 1 each use different methods to solve this problem. The assembly methods of each design should 
be considered to develop the next concept iteration. 

6. Additive manufacturing can reduce the number of cage parts making assembly more 
straightforward. The direct metal laser sintering process used for Con 10b, however, had the 
unintended consequence of warping thin sections of the cage. Alternative additive processes 
should be tried. 

7. Coni lb demonstrated that the compression cage may be created using traditional manufacturing 
by using several individual parts. The required interference fit between the links and hoop walls 
make the structure difficult to assemble. Future multi-part designs should use a method to enable 
clearance between parts during assembly. 

4.3 Plans for Phase 2 

In phase 2, the vibration ring will be optimized for a specific application by the following process. 

1. Determine target vibration ring properties (e.g., stiffness, loss factor, and frequency band) for a 
specific driveline application. 

2. Find the properties of the existing vibration ring prototypes using the load frame test, in order to 
establish a baseline for improved performance and provide model validation data. 

3. Update the analytical framework to include the effect of link deformation, and then validate the 
model using the load test data. 

4. Using the model, define a vibration ring configuration (e.g., electromechanical coupling factor, 
circuit impedance, stack and cage stiffness, and link angle) that would have the desired 
properties. 

5. Iteratively design, build, and test new concepts, based on the desired configuration, using the 
load frame test to evaluate properties. 

6. Evaluate the most promising designs in the driveline application to determine damping 
performance. 
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